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Abstract 

The problem of maximizing weighted sum rates in the downlink of a multicell environment is 
of considerable interest. Unfortunately, this problem is known to be NP-hard. For the case of multi- 
antenna base stations and single antenna mobile terminals, we devise a low complexity, fast and provably 
convergent algorithm that locally optimizes the weighted sum rate in the downlink of the system. In 
particular, we derive an iterative second-order cone program formulation of the weighted sum rate 
maximization problem. The algorithm converges to a local optimum within a few iterations. Superior 
performance of the proposed approach is established by numerically comparing it to other known 
solutions. 

Index Terms 

Weighted sum rate maximization, multicell downlink, convex approximation, beamforming. 

I. Introduction 

For multiple-input multiple-output (MIMO) broadcast channels, dirty paper coding (DPC) is known to 
be the capacity-achieving scheme [1]. However, DPC is a nonlinear interference cancellation technique 

Copyright (c) 2012 IEEE. Personal use of this material is permitted. However, permission to use this material for any other 
purposes must be obtained from the IEEE by sending a request to pubs-permissions@ieee.org. 

This research was supported by Tekes (the Finnish Funding Agency for Technology and Innovation), Nokia Siemens Networks, 
Renesas Mobile Europe, Elektrobit, Xilinx, and Academy of Finland. 

The authors are with the Dept. Communications Eng. and Centre for Wireless Communications, University of Oulu, Finland. 
Email: (ltran, mhanif, atolli, markku.juntti}@ee.oulu.fi. 



2 



and thus requires high complexity. Hence, linear precoding techniques are of practical interest. Herein, 
we consider the problem of weighted sum rate maximization (WSRM) with linear transmit precoding 
for multicell multiple-input single-output (MISO) downlink. Unfortunately, the WSRM problem, even 
for single-antenna receivers as considered in this letter, has been shown to be NP-hard in (21. Although 
optimal beamformers can be obtained using the methods presented, for instance, in O-O, they may 
not be practically useful since the complexity of finding optimal designs grows exponentially with the 
problem size. Hence, the need of computationally conducive suboptimal solutions to the WSRM problem 
still remains. 

Since the WSRM problem is nonconvex and NP-hard, there exists a class of beamformer designs which 
are based on achieving the necessary optimal conditions of the WSRM problem. In fact, this philosophy 
has been used, e.g., in |6]]-|[9]]. Interestingly, in (3j , the authors have numerically shown that the suboptimal 
designs that achieve the necessary optimal conditions of the WSRM problem perform very close to the 
optimal design. In |6), the iterative coordinated beamforming algorithm was proposed by manipulating 
the Karush-Kuhn-Tucker (KKT) equations. However, this algorithm is not provably convergent. In (SJ, 
J9), the WSRM problem with joint transceiver design is solved using alternating optimization between 
transmit and receive beamforming. As we show by numerical results, these methods have a slower 
convergence rate compared to our proposed design. 

In this letter, we propose a fast converging algorithm that locally solves the problem of WSRM for 
multicell MISO downlink. The idea of our iterative beamformer design is based on the framework of 
successive convex approximation (SCA) presented in |[T0l . The numerical results show that the proposed 
algorithm converges within a few iterations to a locally optimal point of the WSRM problem. The 
general concept of the SCA method is as follows. In each step of an iterative procedure, we approximate 
the original nonconvex problem by an efficiently solvable convex program and then update the variables 
involved until convergence. We note that in the context of transmit linear precoding for multicell downlink, 
the SCA method has been used, for example, in (7). Basically, this method is based on convex relaxations 
of the rate function and generally arrives at more complex formulations. By proper transformations, we 
approximate the WSRM problem as a second-order cone program (SOCP) in each step of the SCA 
method. Our numerical results show that the proposed algorithm generally performs better than the 
known approaches, in particular, in terms of convergence rate. 

Notation: We use standard notations in this letter. Bold lower and upper case letters represent vectors 
and matrices, respectively; (.) T represents the transpose operator. C axb represents the space of complex 
matrices of dimensions given as superscripts; |c| represents the absolute value of a complex number. 
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Finally, ||.||2 represents the I2 norm. 

II. Problem Formulation 

Consider a system of B coordinated BSs of N transmit antennas each and K single-antenna receivers. 
The set of all K users is denoted by U = {1, 2 . . . ,K}. We assume that data for the fcth user is 
transmitted only from one BS, which is denoted by G B, where B = {1,2,..., B} is the set of all 
BSs. The set of all users served by BS b is denoted by U^. Under flat fading channel conditions, the 
signal received by the /cth user is 

K 

Vk = h bk , k w k d k + ^2 h^^Wjdj + n k (1) 

i=l,i^k 

where j. G C 1 is the channel (row) vector from BS b, L to user k, w k G C Nxl is the beam- 
forming vector (beamformer) from BS b k to user k, d k is the normalized complex data symbol, and 
n k ~ CAf(0, a 2 ) is complex circularly symmetric zero mean Gaussian noise with variance a 2 . The term 

Yli=i i^k hfti.feWidj in (OQ) includes both intra- and inter-cell interference. The total power transmitted by 

.. 1 1 2 

BS b is Sfce£4 || w fc||2- ^ ne SINR ^ k of user k is 

Ih w I 2 

Ik = L ^ f 72- (2) 

In this letter, we are interested in the problem of WSRM under per-BS power constraints^, which is 
formulated as 

maximize Efc=i «fc log 2 (l + 7fc) 

subject to Y.k&j b \\™k\\l < p b, V6 G B 
where a k s are positive weighting factors which are typically introduced to maintain a certain degree of 
fairness among users. As mentioned earlier, since problem © is NP-hard, the globally optimal design 
mainly plays as a theoretical benchmark rather than a practical solution 01 . Herein, we propose a low- 
complexity algorithm that solves ([3]) locally, i.e, satisfies the necessary optimal conditions of ©. 

III. Proposed Low-complexity Beamformer Design 

To arrive at a tractable solution, we note that following monotonicity of logarithmic function, (f3]) is 
equivalent to 

maximize Y\ k (l + r ) k ) ak 

(4) 



subject to Y,k&A„ Il w fell2 ^ p bi Vfe G £ 



It is straightforward to extend the proposed algorithm to handle per-antenna power constraints at each BS. 



4 



which can be equivalently recast as 

maximize JL tk (5a) 

■Wk,t k 

subject to 7fc > tl /ak - 1, Vfc € W (5b) 
Efcew* llwfclH < A, (5c) 

The equivalence of © and ([5]) can be easily recognized by noting the fact that all constraints in (1551 
are active at the optimum. Otherwise, we can obtain a strictly larger objective by increasing tk without 
violating the constraints. Next, by introducing additional slack variables /3k, we can reformulate © as 

maximize Ilfc^' (^ a ) 

w fc ,t fc ,/3 fc 



subject to h fefcifc w fc > x/t 1 ^ - l&.Vfc G W, (6b) 



Im(h 6fe;fe w fc ) = 0, V/c e W, (6c) 

' \ 1/2 

^ 2 + E^ fe IK,fcWi| 2 J <p k ,VkeU, (6d) 

Efc G wJ|w fc ||i<A, V6Gi3. (6e) 

The equivalence between ([5]) and © is justified as follows. First, we note that forcing the imaginary 
part of hf, k ^Wfe to zero in (l6cT ) does not affect the optimality of © since a phase rotation on will 
result in the same objective while satisfying all constraints. Second, we can show that all the constraints 
in (l6dl ) hold with equality at the optimum. Suppose, to the contrary, the constraint for some k in (l6dl ) is 



inactive. Let us define h — fik/v and t/. = {i] 2 (t 1 f / ak — 1) + where rj is a positive scaling factor. 
Since the constraint ( f6dl ) is inactive, we can choose rj > 1 such that the constraints in (1651 and (|6dl ) 
are still met if we replace by However, such a substitution results in a strictly larger 

objective because > tk for 77 > 1. This contradicts the fact that we have obtained an optimal solution. 
As a step toward a low-complexity solution to the WSRM problem, we rewrite the constraint (l6bl) as 



h-6 fci fcw fc > y/x^Pk, Vk^U (7a) 
x k + l>tl /ak , Vk£U. (7b) 

Again, we can easily see that by replacing d6bl with dTal ) and (TTbl , we obtain an equivalent formulation 
of ©. The reason of doing so becomes clear shortly. Let us define f{xk,Pk) = \fx~k/3k for Xk,f3k > 
and focus on the constraint dTal ) first. Note that f(xk,(3k) is nonconvex on the defined domain, and thus 
dTal ) is not a convex constraint. To deal with nonconvex constraints, we invoke a result of |[T0l which 
shows that if we replace f(xk,Pk) by i ts convex upper bound and iteratively solve the resulting problem 
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by judiciously updating the variables until convergence, we can obtain a KKT point of ([6]). To this end, 
for a given fa for all k, we define the function ifTOl 

G(x k ,p k ,fa)±^Pl + ^- k x k (8) 

which arises from the inequality of arithmetic and geometric means of fafi\ and fa l x k . It is easy to 
check that G(x k , fa, fa) is a convex overestimate of f(x k ,/3 k ) for a fixed fa > 0, i.e., G(x k , f3 k ,fa) > 
f(x k ,fa) for all 4>k > 0. Moreover, when 4> k = ^w^-, it is plain to observe 

f(x k ,fa) = G(x k ,fa,fa) (9a) 
Vf(x k ,fa) = VG(x k ,fa,fa) (9b) 

where V/ is the gradient of /. Obviously if f(x k ,(3 k ) is replaced by G(x k ,/3 k , fa), (P7al ) can be formulated 
as a second-order cone (SOC) constraint as we shown in (II ldl ). 

Now we turn our attention to (f7bl . Recall that trj ak is convex if < a k < 1 and concave if a k > 1, 
and that the optimal solution of (0 stays the same if we multiply all a k 's by the same positive constant. 
Thus, we can force (T7bl to be convex by scaling down a k s in (fj) such that < a k < 1 for all k. 
However, in this case, the constraint x k + 1 > t 1 J ah cannot be directly written as an SOC constraint 

n 

for a k G M++o As our goal is to arrive at an SOCP, we instead scale a k 's in (f3]) such that a k > 1 

for all k and thus t]/ ak becomes concave. Again, in the light of ifTOl . we replace the right side of the 

inequality in (l7bl by its upper bound, which now can be obtained by the first order approximation due 

to the concavity of t]/ ak . Precisely, we have 

J_ J_ J 1 

l k ^ h + -^h - % ) (iu) 

(n) 

where t k denotes the value of variable t k in the nth iteration (i.e., the iteration corresponding to 
Algorithm \T\ described later). In fact, we have linearized t\/ ak around the operating point ij! . With 
(TT0l ). (T7bl now becomes a linear inequality. We note that the linear approximation in (fTOb is trivially 
shown to satisfy the conditions in (|9at and (|9bl at A question naturally arises is whether the linear 
approximation in ( fTOb affects the optimal sum rate. Interestingly, our numerical experiments show that 
the WSR obtained with the successive approximation with a k > 1 is identical to that when (f7bl is forced 
to be convex by having < a k < 1 in (3) for all k. In terms of complexity, the linear inequality in 
(fTOb is more preferable since it requires lower computational effort compared to the original nonlinear 
equality constraint in d7bl . 

2 When a k is an integer or a rational number, we can transform the constraint (T7b) into a number of SOC constraints 1111 . 
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maximize z^ (11a) 

Wk,tk,Xk,/3k,zi 

subject to || [2zl N ~ 1] (t 2 i-i-t 2i )] J \\ 2 < (t 2 i-i+t 2 i), i = l,...,2 N ~ 1 (lib) 



(z?-zW)]\<(4 1) + 4 1) ), (He) 



|| [i(h fefcjfc w fc - ^jx k - 1) V V^] || a < - ^jx k + l),V& G U (lid) 

+ " 4 n) ) < ^ + 1, Vfc € U, (He) 

|| [a h 6l>fc wi • • • h 6fc _ lifc w fc _i h bk+uk w k+ i ■ ■ ■ h bii . ifc w^] T || 2 < f3 k , Vk G U, (1 If) 
E feeW Jw fe ||!<A, V6GB. (llg) 



Replacing the right sides of (TTal ) and (ITbl by the upper bounds in ([8]) and (fTOl ). respectively, we can 
formulate © as an SOCP by noting that the objective in (O, i.e., the product of t k s admits an SOC 
representation ifTTTl . lfT2l . The main ingredient in arriving at the SOCP representation is the fact that the 
hyperbolic constraint uv > z 2 where u > 0, v > is equivalent to \\[2z (u — u)] T ||2 < (u + v). Let us 
illustrate the SOCP formulation of ([6]) for the special case K = 2 q , where q is some positive integer. By 
collecting two variables at a time and incorporating the additional hyperbolic constraint corresponding 
to them, we rewrite © as the SOCP in ([111 ). shown on the top of the page, where is the value of 
<j) k in the nth iteration. In the case of K ^ 2 q , we define additional tj = 1 for j = K + 1, . . . , 2^ og 2 K ~\ , 
where \x\ is the smallest integer not less than x and the above expression still holds iTTTTl . Now we are 
in a position to present an algorithm that solves problem © locally. The pseudocode of the beamformer 
design is outlined in Algorithm Q] We now present the convergence analysis of Algorithm Q] Consider 



Algorithm 1 Proposed beamformer design for the WSRM problem in multicell MISO downlink. 
Initialization: n = 0, {4> k , ijj, n ') = random. 
1: repeat 

2: Solve ([□) with 4>[ n) 

3: Update x^) = {t%,^x%) and ^ 



4: until convergence 




Fig. 1. Average sum rate comparison for single-cell scenario, N = 4, K = 4. 



the n + 1st iteration of Algorithm Q] that solves the optimization problem (fTTI) . If we replace (tk, /3k, Xk) 
by (4 n) )4 n) :4 n) ) and w ^ b y w Mc> a11 tne constraints in flTH)-(lTg) are still satisfied. That is to 
say, the optimal solution of the nth iteration is a feasible point of the problem in the n + 1st iteration. 
Thus, the objective obtained in the n + 1st iteration is larger than or equal to that in the nth iteration. In 
other words, Algorithm Q] generates a nondecreasing sequence of objective values. Moreover, the problem 
is bounded above due to the power constraints. Hence, Algorithm Q] converges to some local optimum 
solution of (fTTI ). By the two properties shown in (© and based on the arguments presented in |[T0l , it 
can be shown that this solution also satisfies the KKT conditions of ©. Numerical results in Section JV] 
confirm that Algorithm \T\ performs very close to optimal linear design. 

IV. Numerical Results 

In this section, we numerically evaluate the performance of Algorithm Q] under different setups using 
YALMIP |[T3l with SDPT3 Ifl4l as internal solver. In the first experiment, we consider a single-cell 
scenario where a BS with N = 4 transmit antennas serves K = 4 users. The entries of h{, & are CAA(0, 1) 
and the noise variance a 2 = 1. In Fig. [T] we plot the average sum rate (a^ = 1 for all k) versus the total 
transmit power P at the BS. The achieved sum rate of Algorithm Q] is compared to those of zero-forcing 
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beamforming lfl~5l . the weighted sum mean-square eiTor minimization (WMMSE) algorithm in the 
soft inference nulling (SIN) scheme in Q, and the optimal linear design using the branch-and-bound (BB) 
method in (3j , 0]. Initial values for beamformers for the suboptimal schemes in Q, 13, and (<^[°\ ij?) 
in Algorithm Q] are generated randomly. The sum rate is obtained after Algorithm Q] and the iterative 
suboptimal schemes in (7], O converge, i.e., the increase in the objective value between two consecutive 
iterations is less then 10 -2 . The gap tolerance between the upper and lower bounds for the BB method 
is set to 10 -1 as in 0, (H, and the resulting sum rate is calculated as the average of the upper and 
lower bounds^ Results reveal that the average sum rate of Algorithm Q] and other iterative beamformer 
designs is the same on convergence and close to that of the optimal linear approach. However, the SIN 
scheme, WMMSE algorithm and the optimal design have a slower convergence rate as discussed next. 

In the second experiment, we illustrate the convergence rate of all considered iterative suboptimal 
schemes. A simple two-cell scenario with each BS serving 2 users is considered. The number of 
transmit antennas at each BS is set to N = 8. The weights, without loss of generality, are taken as 
(ai,a 2 ,a3,a 4 ) = (0.14,0.21,0.28,0.36) and the power budget of each BS is set to P b = 12 dB for 
b = 1,2. Fig. |2] compares the weighted sum rate of the considered schemes as a function of iterations 
needed to obtain a stabilized output for a random channel realization. In particular, our algorithm has 
converged just after a few of iterations, while the WMMSE algorithm is still less than midway to 
convergence. In fact, for this particular case the WMMSE took hundreds of runs before converging 
to the local optimum solution. This observation may be attributed to the fact that optimization strategy 
of requires alternate updates between transmit and receive beamformers and therefore exhibits slower 
convergence properties. We have also noticed that for certain initial values the convergence rate of the 
WMMSE algorithm is greatly improved. In fact, it was reported in |[T6l that the convergence rate of an 
alternating optimization algorithm depends on the initial values of the variables involved, and it converges 
quickly if initial guess is relatively close to the optimal solution. Further, we observe that Algorithm Q] 
with and without scaling has slightly different convergence rate (labeled in Fig. |2]as 'approximated' and 
'not approximated') but same optimal value. This validates that the approximation used to arrive at an 
SOCP formulation has no impact on the achieved sum rate. Our numerical results reveal that for other 
channel realizations, the SIN scheme may have similar convergence behavior to Algorithm [T] but the 
average per iteration running time of Algorithm Q] is approximately four times less than that of the SIN 

3 The principle of the BB method to compute an optimal solution to nonconvex problems is to find provable lower and upper 
bounds on the globally optimal value and guarantee that the bounds converge as iterations evolve. 
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1 3 5 7 9 11 13 15 17 19 

Iteration index 

Fig. 2. Convergence rate of the weighted sum rate, B — 2, N = 8, K = 4, Pi = 12 dB for all b. 

method. For the set of channel realizations considered in Fig. |2j the optimal design also converges to the 
same point achieved by other iterative suboptimal methods. However, it takes more than 600 iterations to 
reduce the gap between lower and upper bounds of the BB method to less than 10 _1 . Theoretically, the 
faster convergence of Algorithm Q] may be attributed to solving an explicit SOCP in each of its iterations. 
The faster convergence of our algorithm can be much useful for distributed implementation which is left 
as future work. 

V. Conclusion 

In the letter we have studied the problem of WSRM in the downlink of multicell MISO system. 
Since the problem is NP-hard, we have proposed a low-complexity approximation of the optimization 
problem. We show that the problem can be approximated by an iterative SOCP procedure. While the 
convergence of the algorithm can be proved, its global optimality cannot be established. Nonetheless, the 
algorithm outperforms the previously studied solutions to the WSRM problem, in particular, in terms of 
its convergence rate. 
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